o 



X 



Relativistic hydrodynamics with strangeness 
production 



L Turko 



o 

o . 

, Institute of Theoretical Physics, University of Wroclaw, 

^ I pi. Maksa Borna 9, 50-204 Wroclaw, Poland 



E-mail: turko@ift.uiii.wroc.pl 



' Abstract. The relativistic hydrodynamic approach is used to describe production of 

strangeness and/or heavy quarks in ultrarelativistic heavy ion reactions. Production 
f— I ^ processes are important ingredients of dissipative effects in the hadronic liquid. Beyond 

Q^' viscosity also chemo- and thermo-diffusion processes are considered. This also allows 



to specify chemical and thermal freeze-out conditions. 
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^ ' 1. Introduction 

'NT 

^ . Results of RHIC experiments showed up that data of ultrarelativistic heavy ion collisions 
. can be described using relativistic hydrodynamic of a thermalized perfect fluid [H [2]. 
Indications came mainly from a large elliptic flow effect. It appears, however, that 



dissipative effects become more and more important [3j beyond mid-rapidity, most 
central, and the top RHIC energy regions. A more detailed analysis of initial conditions 
of color glass condensate models needs a finite viscosity to reduce the elliptic flow to 
match the data |4j. There are also some general arguments that at least shear viscosity 
of matter created in heavy ion collision must be nonzero [3, E] . 

The entropy of the QGP phase at hadronization determines the observable particle 
multiplicity [3, [H |9] so it is important to take into account contributions due to 
particle production processes. The entropy content is usually estimated assuming 
thermal momentum distributions and related equations of state. A more general 
formalism allowing for the consistent incorporation of dissipative and entropy producing 
mechanisms into the relativistic hydrodynamical framework was presented in |10j . 
That was the hydrodynamic model with dynamical chemical reactions and multi-fluid 
included. We are going to specify here this general scheme for the case of strangeness 
production in the quark-gluon plasma (QGP) case. 

The entropy balance is treated here in the deconfined QGP phase. The primary 
contribution is from the initial decoherence processes which however cannot be 
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calculated ab initio. This contribution is estimated from 40 [TT] to 70 % [12] of the total 
entropy produced. We consider the contributions due to the production of secondary 
partons. This is the main calculable entropy source. The entropy originated from the 
hadronization processes is rather small. 

The energy-momentum tensor for an imperfect fluid is obtained in the standard 
procedure [13]. The energy- momentum tensor is given as 

T^"" = (e + P)m^m'^ - Pg^"" + ST^"" , (1) 

with the term 6T^'^ related to all dissipative processes in the fluid, e and P denote the 
local energy density and pressure, respectively. 

The evolution of the closed system is described by equations 

d^T"" = . (2) 

These equations are supplemented by continuity equations 

= d^p.u^ = , (3) 

where pi denotes the particle number density in the local rest frame. The local rest 
frame is determined by the four-velocity condition = (1, 0, 0, 0). 
We require dissipative terms to vanish in the local rest frame, 

Tj"] = e 5T°° = ^ u,u,5T'^'' = . (4) 

These conditions allow to determine the general structure of the dissipative term 6T'^'^ 

m. 



1.1. Diffusion, particle production, chemical reactions 

We generalize continuity equations ([3]) to incorporate diffusion and particle production 
or chemical reactions. Equation ([3]) is replaced now by 

d,pf = + A'^^a,7^,) = J, . (5) 

with Af^i^ = g^u — Uf^Uu, g^^ = diag(l,-l,-l,-l) the metric of the fiat space-time. 

TZ is the diffusion term which is determined consistently [10] with the Second Law 
of Thermodynamics. J7i are sources of particle production. They are determined by 
microscopic reaction cross sections. In the local rest frame one gets the gain and loss 
equation 

j 

Since the entropy of the QGP phase at hadronization determines the observable 
particle multiplicity, it is crucial to understand its production. The entropy production 
in a dissipative relativistic fluid is given by: 

i 

+ Y^ip, + TOd^A^^d^TZi + T^(9^A)A'^^a,7^, , (7) 
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where T is the temperature, = 1, and we employed natural units h = c = ks = I- 
The heat-flow four- vector is defined by 

Q^ = d^T-Tu''d,u^, (8) 

and the shear tensor is given by 

2 

W^i. = d^u^ + duU^ - -g^^d^u'^ . (9) 

2. Entropy Production 

The First Law of Thermodynamics imply: 

= Ta^S^ + J2 - d^A'^'^dM + u^d^ST'^" , (10) 

i 

with an auxiliary entropy four-current, s'^ = su^, and = piu'^ 

Chemical potentials /ij do not imply here traditional chemical equilibrium as they 
are space and time dependent parameters here. They are related to abelian charges 
conservation such as the electric charge, the strangeness or the baryon number. In the 
case of a conserved abelian charge the corresponding chemical potential is written as: 

i 

where denotes the charge of particles of species " i" . 

Taking into account of nonabelian symmetries conservation is more involved and 
can be formulated on the kinetic equation level [TH [T^ . 

The proper entropy four-current is 

s^" = s^ + T-\^6T^"' + Y,^iA^''di^T^i ■ (H) 

i 

2.1. Second Law of Thermodynamics 

An increasing entropy conditions leads to constraints on the particle source terms J7i 
and on the diffusion terms 71 and £ [lOj . 

The system is driven towards chemical equilibrium if 

Y.liiJi<Q . (12a) 

i 

This condition is closely related to abelian charge conservation of the system. 

The structure of the dissipative term 5T^^ compatible with d^j^s^ > is given as 

5T^''' = /€(A^^m'^ + A'^^M^)Q^ + r]A''^'A''^W^s + (A'^'d^u^ , {12b) 

where k, rj and ( denotes the coefficients of heat conductivity, shear viscosity, and bulk 
viscosity, respectively. The heat-fiow and the shear tensor terms are given here by ([8]) 
and (ED. 
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Also, the contribution to the entropy production involving TZ becomes a nonnegative 
if one identifies: 

A = , 7^.^E^^.f , (12c) 

i 

The symmetric matrix a consists of the mutual diffusion constants and has only 
nonnegative eigenvalues. 
The diffusion current 

I ^ -V7^, = - E ^(V/i, - |VT) . (13) 

i 

involves chemo- and thermo-diffusion contributions. 

The entropy production in the relativistic fluid is now given as 

- ^Q.A^'-Q. + |wf A^W^ A? + C{d,un' . (14) 

This allows us to specify conditions for particular freeze-outs encountered in high 
energy heavy ion collisions. The chemical freeze-out is reached with the hadrochemical 
equilibration. An universal chemical equilibrium means that all sources JT^ = 0. A 
partial chemical equilibrium is also possible when only some of sources vanish. That 
would correspond to the situation when e.g. light quarks equilibrate faster than heavy 
quarks. We have then different chemical freeze-outs. 

A complete thermal equilibrium is reached when the entropy is not longer produced. 
It just means d^s'^ = condition. 

3. Strangeness Production 

We are going to specify terms related to the strangeness production in heavy ion collision 
processes. Strangeness production in plasma is dominated by gluon fusion g + g s + s, 
with much smaller contributions of light quarks annihilation processes q + q —>■ s + s. 
We specify the gain and loss equation ([6]) in the local rest frame 

dtPi ^J^ = Gs-Ls. (15) 

The gain term is given as 

Gs = {crl^M' + E {^!%,)p,p, ^ g . (16a) 

q=u,d 

The loss term is 

The reactions rates {cTabVab) are cross sections averaged over momentum distributions of 
the components for the gain and the loss process, respectively. The appearance of the 
pI term here instead of the psPs term reflects the local strangeness conservation, ps is 
now the density of produced ss pairs. 
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The source term for strangeness has now a standard form \W\ of the kinetic rate 
equation 

^ = g-cpi (17) 

Gluons and hght quarks numbers are considered to be abundant so we can neglect 
the change of their numbers due to the strangeness related processes. It means that 
gluons and hght quarks are, in a sense, treated as a heat bath for strange quarks. Similar 
situation appears also for other heavy quarks processes. 

It appears, however, that local exact strangeness conservation leads to strong 
correlations which should be taken into account [17j when Equation (fT7|) is analyzed. 
The density ps in this equation is in fact an event-by-event average density. Let's denote 
event-by-event averages by ((■)). The rate equation f|T71) is changed to 

^iM = g_ cUpD) =g- C{{p,))' - C6pl , (18) 

where is the event-by-event fluctuation of the number of ss pairs. 

The relation of this fluctuation to the event-by-event average of the number of pairs 
density if crucial for the behavior of the ss pairs density. If ((ps))^ ^ Spl then 

Ps(r) = ^tanh(rv/^) , 

- a characteristic behavior for the grand canonical distribution. 

Increasing fluctuations are related to the exact local strangeness conservation. This 
effect becomes more profound for strangeness-poor systems. The system would behave 

- in the simplest scenario [TT] - according to the canonical distribution. The relation 
canonical - grand canonical description is more involved in practical applications. It 
appears [TS] also that geometric effects should be taken into account as the strangeness 
distribution is centrality dependent. 



4. Conclusions 



We considered here the entropy producing mechanisms of chemo- and thermo-diffusion 
together with the contribution of particle production. It was reached within the first 
order formulation of diffusion process which usually leads to some difficulties. The most 
serious one is the instability of homogenous equilibrium when even small perturbations of 
the heat flow, flux velocity or the internal energy density lead system out of equilibrium. 
It was pointed out recently [I9] for perfect fluid that by a particular separation of 
internal and flow energies such instabilities can be removed. It seems also that viscosity 
and particle production processes added to the standard diffusion processes stabilize a 
theory and improve the situation even further. 

The presented model of strangeness production can be directly extended on other 
heavy quarks processes. It would be interesting to estimate equilibration time-scales for 
different flavors at different energies. An implementation of the appropriate numerical 
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code with all initial conditions subtleties is a separate topic. This is necessary to get some 
quantitative estimations of considered effects. Recent numerical results [20] of the 2+1 
shear-viscous one-component hydrodynamics show the quantitative importance of the 
viscosity in the analysis of the system. Cooling is slower and the transverse expansion 
is faster. It should be stressed that an inclusion of production processes makes the 
problem more involved. As was shown |12] quantitatively, within the thermal model, 
strangeness production leads to the faster coohng of the fluid. This is related to the 
conversion of the initial condition pressure into strangeness. This effect should be taken 
account within the hydrodynamical model. The resulting "effective viscosity" decreases 
so the fluid becomes more and more "effectively perfect". 

An interesting point is the possible dynamical behavior of bulk viscosity terms. The 
relaxation time for strange and heavy quarks is substantially larger than for light quarks. 
This increases equilibrium- nonequilibrium "amplitude", resulting in the increasing bulk 
viscosity term. 

These problems are objectives of our future research plans. 
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